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CHAPTER  I 
INTRODUCTION 

If  one  accepts  the  unif ormitarian  theory  of  the 
evolution  of  the  earth,  then  earthquakes  have  occurred 
throughout  the  long  history  of  the  earth,  as  they  do  at  the 
present  time.  The  recorded  evidence  for  the  occurrence  of 
earthquakes  only  dates  back  to  a  few  millennia.  Branson  et  al % 
(1955) ,  make  the  mention  of  earthquake  occurrences  as  reported 
by  historians  such  as  Pliny  and  Herodotus  dating  back  2500  years. 
That  earthquakes  release  energy  and  give  rise  to  elastic  waves 
has  been  known  only  for  the  last  few  centuries  and  detailed 
studies  of  the  propagation  of  these  elastic  waves  have  only 
been  carried  out  relatively  recently,  having  been  begun  by  such 
workers  as  Rayleigh,  Love,  Lamb,  and  Jeffreys.  Seismic  waves 
generated  by  earthquakes  have  been  used  extensively  to  study  the 
details  of  the  crust  of  the  earth.  Mohorovi£id  was  one  of  the 
first  to  demonstrate  the  existence  of  a  crust  of  lower  velocity 
material  overlying  a  higher  velocity  mantle. 

One  of  the  major  difficulties  in  utilizing  earthquake 
waves  to  study  crustal  structure  is  the  uncertainty  in  both 
time  and  location  of  the  energy  release.  Thus  in  recent  years 
investigators  have  come  to  rely  more  and  more  on  man-made 
explosions  where  both  the  location  and  time  of  the  explosion 
are  accurately  known. 
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Geologists  have  shown  that  the  sedimentary  rocks  near  the 
surface  of  the  earth  are  often  layered  with  a  particular  rock  unit  a 
few  tens  of  feet  thick  often  occurring  over  an  area  of  many  tens  or 
even  hundreds  of  square  miles.  Thus  the  concept  of  plane  horizontal 
layers  is  a  justifiable  one  near  the  surface  of  the  earth.  Travel¬ 
time  curves  for  the  acoustic  waves  from  earthquakes  and  man-made  explo¬ 
sions  often  suggest  that  the  earth  is  layered  even  below  this  upper 
sedimentary  section  and  thus  the  assumption  of  layering  of  the  crust  as 
a  whole  is  a  common,  and  useful  approximation. 

In  an  elastic  solid,  two  types  of  waves  may  be  propagated,  lon¬ 
gitudinal  or  P-wave,  and  shear  or  S-wave.  There  is  a  considerable 
body  of  work  (an  extensive  list  of  references  is  given  by  Ewing  et  al .  , 
1957)9  both  theoretical  and  experimental,  dealing  with  the  reflection 
and  refraction  of  these  elastic  waves  at  the  boundaries  of  various 
layers.  Waves  incident  on  an  interface  with  the  source  of  disturbance, 
say  in  the  first  medium,  are  refracted  and  reflected  at  the  interface 
depending  on  the  boundary  conditions.  In  the  case  of  acoustic  waves 
one  often  makes  use  of  the  critical  angle  as  used  in  optics.  When  the 
angle  of  the  incident  wave  is  equal  to  the  critical  angle,  the  elastic 
waves  travel  along  the  boundary,  and  are  then  returned  into  the  medium 
in  which  the  source  is  located,  as  so-called  head  waves.  These  waves 
can  undergo  conversion  at  the  interface  from  P  to  S  or  vice  versa  de¬ 
pending  again  on  boundary  conditions.  The  present  investigations  deal 
with  these  converted  phases. 


3 


The  crustal  structure  of  Southern  Alberta,  based  on  the  refrac¬ 
tion  waves  generated  by  explosions 9  has  been  the  subject  of  several 
investigations  (Weaver  1962,  and  Maureau  1964 )9  and  these  results  have 
been  reported  in  papers  by  members  of  the  Physics  Department,  University 
of  Alberta  (1962,  1964,  1966).  Fig.  I.  1  shows  the  location  of  various 
profiles  shot  in  the  general  area  of  Southern  Alberta  and  Saskatchewan. 
Bow  City  Shot  Point  and  Suffield  Shot  Point  locations  can  be  seen  in 
this  map. 

The  present  work  is  based  on  data  from  two  refraction  profiles 
in  Southern  Alberta.  One  profile  utilizes  recordings  of  seismic  waves 
produced  by  a  500  ton  charge  of  TNT  exploded  on  the  surface  at  the 
Suffield  experimental  station  of  the  Defence  Research  Board,  and  re¬ 
corded  in  the  foothills  and  Rocky  Mountains;  the  other  consists  of 
recordings  from  a  series  of  buried  charges  of  various  sizes  up  to  1000 
pounds,  at  a  shot  point  near  the  town  of  Bow  City.  The  detector  loca¬ 
tions  for  the  second  profile  are  in  approximately  the  same  place  as 
those  for  the  first  and  thus  direct  comparison  can  be  made  between  the 
data  on  each  profile.  Unfortunately  the  profiles  are  shot  in  the  same 
direction  and  this  results  in  an  uncertainty  in  the  dip  of  the  various 
layers  from  which  arrivals  have  been  detected.  However  the  records 
examined  furnish  excellent  opportunities  for  the  study  of  later  events 
which  throw  considerable  light  on  the  structure.  In  principle  it 
should  be  possible  to  resolve  the  uncertainties  of  dip  by  utilizing 
converted  phases,  but  in  practice  this  appears  to  be  difficult. 
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CHAPTER  II 


CONVERTED  PHASES 


II.  1  General  Considerations 

II.  1.  1  The  two  types  of  elastic  waves,  compressional  or  P,  and  shear 
or  S,  are  known  as  body  waves.  An  explosion  or  an  earthquake  releases 
energy,  some  of  which  is  propagated  in  the  rigid  earth  by  these  waves. 
When  these  waves  are  incident  on  an  interface  they  are  either  refracted, 
or  reflected  at  the  interface.  For  two  layers  in  "welded"  contact  at 
z  =  0,  with  the  source  of  disturbance  in  layer  1,  the  displacements  are 
given  (Ewing  et  al.  1957),  by  the  following  equations. 


for  i  =  1,  2 


(1) 


where  the  potentials  <f)^  and  iJk  are  solutions  of  the  equations 


for  i  =  1,  2 


(2) 
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Ewing  al*  give  the  following  solutions  for  the  above 
equations 


00  k  —  v !  |  z— h  |  oo  -v  (z-h) 

*1  =  /  VT  J0(kr)e  dk  +  /  Q,(k)J0(kr)e  dk  (3) 

0  V1  o  A 


oo  -vj(z-h) 

=  f  S1(k)J0(kr)e  dk  (U) 


(3)  and  (b)  hold  for  z  >  0 

oo  v  (z-h) 

<f>2  =  /  Q2(k)J0(kr)e  dk  (5) 


oo  v2(z-h) 

=  /  S2(k)Jc(kr)e  dk  (6) 

(5)  and  (6)  hold  for  z  <  0 


v;  = 

/?- 

k62! 

(T) 

^2  = 

*  /  k2 

"  k3z 

(8) 

V1  = 

J  k2  - 

k  2 
al 

(9) 

v2  = 

/  k2  - 

k  2 
ka2 

(10) 

The  four  coefficients  Qj_ ,  can  be  chosen  to  satisfy  the  boundary 


conditions . 


. 
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At  z  =  0 , 

<*!  =  q2  and  Vl  =  W2  (ll) 

^zz^i  ~  ^zz^2  i  “  ^zr^2 

(ll)  merely  state  that  the  displacements  and  stresses  are  continuous 

across  the  interface. 

Hence  we  have 


3<j>,  3<J>2  32^2 

9r  +  3z3r  3r  +  3z3r 


(12) 


3<f>i  32^x 

3z  3z2 


+  ko2^i  = 

Pi  1  3z  3z2  P2  z 


3(J>2  3  ^2  2 


(13) 


9  3w  3w 

\v2*l  +  2U1  3^=  X2V ^2  +  2M2  H 


(Ht) 


3q.  3w_  3q  3w0 

**1  (^3^ 


(15) 


Since  <f>  1 ,  t|>. ,  <p2  and  \p2  are  given  by  (3),  (U),  (5)  and  (6)  above,  one 
can  differentiate  these  functions  with  respect  to  z  and  r»  and  by 
using  the  properties  of  Bessel's  Functions  one  can  show  that  (12)  to 
(15)  are  reduced  to 


Ql  -  q2  -  v jS !  -  v2S2  = 


-v  ih 


vi 


(16) 


vl^l  +  V2Q2  “  k2Si  +  k2S2  =  ke 


-vih 


(17) 
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-  (2y1k2  -  p1w2)Q1  +  (2y2k2  -  p2v2)Q2  +  2y1k2vjs1 


k  -vjh 

+  2y2k2v2S2  =  —  (2yxk2  -  p1w2)e 


(18) 


2y i vxQi  +  2y2v2Q2  “  (2yxk2  -  PjW2)^  +  (2y2k2  -  p2w2)§2  =  2kyxe 


— v  j  h 


where 


vih 

Qi  =  Qie  i  Q2  =  Q2e 


-v2h 


(19) 


i ,  « , 

vlh  *  ‘v2h 

Si  —  S^e  j  s2  -  S2e 


If  we  now  put  (Ewing  et  al . ) 


2yxk2  -  PiW2  =  ax  ;  2y2k2  -  p2w2  =  a2  , 

then  the  determinant  of  the  above  systems  of  equations  (l6)  to  (19) 
is  given  by 


A(k)  = 


V1 
-ai 
2y  i  vi 


-l 

v2 

a2 

2y2v2 


-k‘ 


2yxk2vJ 


•8.1 


_v2 

k2 

2y2k2v2 

a2 


Using  Kramer's  Rule  we  can  now  write  down  the  coefficients  and 
as  follows . 

a;  -(vi+vj)h 

7 


*1  =Te 


Ax  “2vlh 


u  i 

Sx  -  — i-  e 
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where 


A1  = 


Q2  =  ~  e 


-1 


A2  “(vi“v2^h 


2y0v 


2  2 


-V 


-k 


1 

2 

.2  ' 


-a. 


S2  ■  — « 


ij  -(Vj-v^Jh 


-V 


2y1k2v/1  2y2k2v"2 


k 

v. 


ka 


1 


1 

2kyx 


A2  " 


-k" 


■v, 


2y1k  vj  2y2k2v2 


-k/vj 


ka! 


2ky . 


2y  i  v 


1  V1 


9  9 

Two  similar  expressions  can  be  given  for  Aj  and  A2„  Therefore  we 
write  down9  following  Ewing  et  al. a  (3)9  (*0»  (5)  and  (6)  as 


*i  =  l~  J0(kr)' 


-Vj  |  z-h  |  oo  a  -v  (z+h) 


o  v 


dk  +  /  T"  uo 


(kr  )e 


dk 


_a  -ViZ-Vih 

=  j  “  JQ ( kr ) e  dk 


A1 


°o  A  V2z-Vih 

<f>2  =  /  ~  J0^kr^e  ^ 

o 


can 


(20) 


(21) 


(22) 
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^2  =  /  T-  Jo<kr)e 


v2z-v1h 


dk 


(23) 


The  first  term  in  (20)  represents  the  compressional  wave.  All  other 
terms  in  (20)  to  (23)  represent,  according  to  Ewing  et  al .  ,  waves 
generated  by  it.  To  investigate  these  waves  one  has  to  insert  the 
time  factor  and  carry  out  the  integration.  Different  types  of  waves 
are  determined  by  a  set  of  branch  line  integrals  corresponding  to 
K  =  Ka^ ,  Ka^  and  and  by  residues  corresponding  to  roots  of 

the  equation 


A  (k)  =  0 


Or  A(k)  can  be  given  by 


w2(p2  -  pj) 


A(k)  =  o  =  My.  -  yj2  {[k2(k2  -  — ■ 

2  1  2(y  -  Ul) 


-)2 


'ivJCk2 


p2w" 


2T 


Wl) 


plw 

"  V 


v2v'[k2  +  ^rrL— -T]2 


(V1V2  +  V2V1}  *  VlV2VlV2k2> 


K  has  four  roots  and  accordingly  we  have  four  cases,  each  case  corre¬ 
sponding  to  a  wave  type. 

One  can  consider  the  wave  path  as  composed  of  three  parts.  One 
from  source  to  interface,  second  along  the  interface  and  third  from 
interface  to  detector.  The  coefficient  of  h  in  the  exponential  states 
whether  the  first  part  is  compressional  or  shear  wave.  The  coeffi¬ 
cient  of  z  in  the  exponential  gives  the  same  information  about  the 
third  part,  while  the  value  of  ’k*  at  the  branch  point  indicates  the 
mode  of  travel  along  the  interface.  If  we  have  two  media  with  com- 


i  at  >  (  jrf)A  ?0 
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pressional  wave  velocities  V1 ,  and  V2  and  shear  wave  velocities  v^ ,  v2 
and  such  that  V2  >  v2  >  V1  >  v^  then  there  is  a  possibility  of  shear 
wave  along  the  interface  as  well.  In  the  case  of  the  earth  one  rarely 
comes  across  two  media  such  that  the  shear  wave  velocity  in  the  lower 
medium  is  greater  than  the  compress ional  wave  velocity  in  the  upper 
medium,  so  that  if  the  wave  travelling  along  the  boundary  of  the  lower 
medium  is  a  shear  wave  it  cannot  give  rise  to,  or  be  induced  by,  a 
compressional  wave  in  the  upper  medium.,  This  limits  the  number  of 
possible  cases  to  the  wave  types  shown  in  Figure  II ,1.  If  the  layer 
boundary  is  horizontal,  2)  and  3)  will  arrive  simultaneously  and  the 
only  means  of  distinguishing  between  them  will  be  an  analysis  of  the 
particle  motion  at  the  detector. 

Waves  of  types  l)  to  U)  will  have  the  same  velocity  across  the 
spread  of  detectors  if  there  is  no  dip,  but  those  waves  having  one  or 
more  paths  in  the  upper  medium  as  an  S  wave,  will  be  delayed  with 
respect  to  the  P  wave  by  an  amount  which  depends  on  Poissons' s  ratio 
in  the  upper  layer  and  on  the  thickness  of  the  upper  layer.  The 
analysis  of  crustal  structure  is  based  on  these  time  differences. 

In  general  there  may  be  several  layers  above  the  refracting 
horizon,  so  there  will  exist  the  possibility  of  further  conversions 
at  each  of  the  interfaces.  Boundary  conditions  will  determine  the 
relative  amounts  of  energy  in  each  of  these  possible  waves  and  hence 
will  determine  whether  any  particular  wave  type  may  be  detected  at 


the  surface. 


(4) 


(5) 
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II.  1.  2  Hall  (1964)  derived  equations  for  the  arrival  times  of  waves 
for  various  cases  of  conversions  at  gently  dipping  interfaces.  In  a 
later  paper  (1966)  the  same  author  derives  the  equations  for  the  con¬ 
verted  phases  at  steeply  dipping  interfaces.  The  case  dealt  here  is 
that  of  gently  dipping  beds.  Hence  all  the  calculations  made  in  the 
present  work  are  based  on  the  derivations  given  in  Hall's  first  paper 
(op.  cito)  Hall  and  Brisbin  (1965)  have  used  the  same  theory  to  explain 
various  events  observed  on  records  from  a  profile  at  Flin  Flon- 
Mafeking,  Manitoba.  Pakiser  and  Hill  (1963)  used  later  events  to 
interpret  the  crustal  structure  in  Nevada  and  Southern  Idaho.  Cook  et 
al,  (1962)  have  discussed  the  status  of  PS  converted  phases  in  crustal 
studies  . 

II.  1.  3  A  PS  converted  wave  (Fig.  II.  2)  is  a  seismic  body  wave  which 
results  from  the  conversion  of  an  incident  P-wave  at  a  boundary  within 
the  crust,  to  an  S-wave.  This  S-wave  must  be  a  vertically  polarized 
shear  or  Sy wave  in  order  to  satisfy  the  boundary  conditions. 

An  SP  converted  wave  is  a  seismic  body  wave  which  results  from 
the  conversion  of  an  incident  S-wave  at  a  boundary  within  the  crust, 
to  a  refracted  P-wave. 

In  order  to  positively  identify  an  Sv  wave  it  is  important  to 
note  that  horizontal  motion  geophones  must  be  employed.  Following 
this  several  other  criteria  may  be  used.  One  of  them  is  the  study  of 
particle  motion  to  confirm  that  these  waves  have  motions  consistent 
with  Sv  waves .  The  second  is  the  study  of  the  amplitude  relationships , 


o 


RECEIVER 


SOURCE 


RECEIVER  SOURCE 


V,  V2  *  LONGITUDINAL  WAVE  VELOCITY 
V,  »  SHEAR-WAVE  VELOCITY 


Fig.  II.  2 


•  .  ' 

YTID0J3V  3 VAW  -  ^  v,3 


-> 


12 


II.  2  Amplitude  Relationships 

II.  2.  1  Theoretical  studies  of  the  refraction  of  elastic  waves  and 
their  conversion  have  been  carried  out  in  several  different  ways.  For 
example  the  relative  partition  of  energy  for  the  case  of  an  incident 
compress ional  wave  or  the  same  for  an  incident  shear  wave,  was  de¬ 
rived  "by  Knott  (1899)  and  the  partition  of  amplitudes  was  expanded  by 
Zoeppritz  (1919). 

Macelwane  and  Sohon  (1932)  give  the  various  expressions  derived 
by  Knott  and  Zoeppritz.  The  reader  is  referred  to  this  excellent  work 
for  details. 

There  are  various  other  works  published  by  different  authors 
such  as  Jeffreys  (1926),  Muskat  (1933),  Muskat  and  Meres  (19^0), 
Gutenberg  (19UU),  Heelan  (1953),  Nafe  (1957),  Schwind  (i960), 

Richards  (1961)  and  Steinhart  et  al.  (1961),  all  of  which  deal  with 
amplitudes  of  incident  P-  or  S-waves  and  the  other  waves  generated 
by  them. 

II.  2.  2  Muskat  and  Meres  (oc.  cit.)  start  with  the  wave  equations 
of  the  various  conversions  at  a  boundary  and  deduce  the  amplitude 
ratios.  The  method  is  briefly  described  below. 

The  wave  equations  are 

(V2  +  h2)$>  =  0  h2  =  w2/a2 

where 

(V2  +  k2)'!'  =  0  k2  =  w2/B2 
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w  being  the  wave  frequency  and  a  and  £  the  longitudinal  and  transverse 
wave  velocities  respectively. 

For  an  incident  longitudinal  wave,  they  give  the  following  system 
of  equations. 


$  =  I  e 

1 


i(a]z-bx)  i(-a.1z-bx) 


+  R^e 


*1  =  1S!e 


i(a2z-bx) 


$  =  I  e 

2  2 


=  iT2e 


i(a2z~bx) 


i(c2z-bx) 


For  an  incident  transverse  wave  the  following  relations  hold. 

i(-a  z-bx) 


*i  =  Rie 


Y  =  iTie 


i(c1z«bx)  i(-c3z~bx) 


+  iSTe 


$  =  I  e 

2  2 


i(a2z-bx) 


i(c9z~bx) 


4'.  =  iT_e 
2  2 


where 


a2  +  b2  =  h2 


c2  +  b2  =  k2 


,  ,  ,  w  sin  0 

b  =  b  =  b  =  - 

12  a1 


In  Fig.  II.  3  are  shown  these  two  cases  when  the  incident  wave 
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is  longitudinal  and  when  it  is  transverse,  and  also  shown  are  the 
various  waves  generated  by  them.  Then  they  go  on  to  derive  the  expres¬ 
sions  for  the  amplitude  ratios  (Sj/Tj),  (i^T^,  and  (T^T^  ) 

where  R^  and  S1  are  reflection  coefficients  and  I2  and  T2  are  transmis¬ 
sion  coefficients.  To  get  the  longitudinal  wave  amplitudes  it  is  only 
necessary  to  multiply  these  coefficients  (i  or  R)  by  h  and  to  get  the 
transverse  wave  amplitudes  it  is  only  necessary  to  multiply  these 
coefficients  (S  or  T)  by  k. 

The  following  relationships  are  readily  derived. 


vr/v 


1 


(ci  -  c2<S) 
(Cj  +  c"2V)_ 


Vv 


1 


2ei 

(ci  +  c?_6) 


From  which  it  can  be  seen  that 


+  1 


where 


vr  ■+  displacement  amplitude  of  reflected  wave 
v^  ->■  displacement  amplitude  of  incident  wave 
v2  displacement  amplitude  of  transmitted  wave 

6  y 2. / P i  being  rigidity  modulii  of  the  two  media. 


Considering  the  wave  intensities  and  energy  partition  it  can  be  easily 


shown  that 
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1  -  By  /(b2  +  (l  -  B2)  sec20) , 2 

1  +  By  /(B2  +  (l  -  B2)  sec2©)' 


and  hence  the  transmission  coefficient  is 


(1 


In  the  above  expression 

5  =  82/0i 

Y  =  Y2/Y1 

Yj»  Y2  are  the  densities  of  two  layers. 

B1S  B2  are  wave  velocities  in  these  layers. 

In  case  of  head  waves  the  angle  becomes  critical  angle  which  can  now  be 
given  as 

a  3  <*2  >  04 

Sin  0  =  Sin  i_  =  — ^  “  for  all 

0  “2  b2  e2  >  8l 

Muskat  and  Meres  (op.  cit.)  have  also  published  numerical  results. 


II.  2,  3  Heelan  (1953)  gives  an  erudite  and  elaborate  treatment  on 
the  head  waves,  and  on  the  incident  P-waves  and  S-waves ,  and  the 
conversions  occurring  thereby. 

Heelan' s  argument  starts  almost  at  the  same  place  where  we  ended 
the  section  1  of  this  chapter.  He  recognizes  a  localized  disturbance 
in  the  upper  medium  radiating  P,  SV  and  SH  waves.  When  these  waves 
are  incident  on  the  interface  several  systems  of  head  waves  and  second 
order  boundary  waves  are  generated,  each  associated  with  grazing 
incidence  of  one  or  the  other  of  the  reflected  or  refracted  waves. 


JSji 
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Associated  with  the  grazing  incidence  of  P1P2>  the  refracted  P-wave, 
is  the  head  wave  system  comprising  P1P2P1 ,  P1P2S1  in  the  upper  medium 
and  P1P2S2  in  the  lower  medium.  Associated  with  the  grazing  incidence 
of  SV-  disturbance  are  Sj^Pi,  S1P2S1  in  the  upper  medium  and  8^282 
in  the  lower  medium.  Either  set  of  the  three  waves,  mentioned  above, 
satisfy  the  conditions  of  continuity  of  displacement  and  stress  across 
the  boundary. 

The  primary  incident  radiation  generated  by  the  source  has  the 
form  of  three  auxiliary  wave  functions 


*o  =  /  e 


00 


ikVt 


(l)  a(z-d) 

dk  /  f0H0  (ar)e  da 


0 


c 


for  longitudinal  incident  wave 


*  ikVt  (1)  6(z-d) 


0Q  =  /  e  dk  /  g0H0  (ar)e  da 


0  c 


for  incident  SV  wave 


and 


00  ikVt 
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(l)  B(z-d) 
dk  /  n0H0  (ar) 


da 


c 


for  incident  SH-wave 


where 


kV  =  hv  and  c  is  a  loop  (°°i,  -k,  -h,  °°i) 


where  arg.  a  =  arg.  a  =  arg.  3  =  tt/2  initially,  and  2t r  ^  arg.  a  ^  0 


. 


••  •  - 


r. 


. 


IT 


on  the  path  , 

fQ  =  P1  (k)Aa(2a2/h2  +  1  -  2v2/V2 ) /87ry  (a2  -  k2)1/2 

g0  =  Pj_  (k)AaATTph2 

n0  =  S1  (k)Aa/Urry (a2  -  h2)1/2 
00 

vhere  P,(t)  =R i  j  P, (k)exp(ikVt )dk 
1  o00  1 

oo 

Sj(t)  =  RJl  j  S1  (k)exp( ikVt )dk 
o 

R£  indicating  the  real  part  of  the  integrals » 

A  =  volume  of  the  cylindrical  source 
y  =  rigidity  of  the  medium 

V,v  =  velocities  of  P  and  S  waves  respectively 0 
Heelan  then  lets  0  and  x  satisfy  the  well  known  wave  equations, 
and  expresses  three  functions  for  both  media  in  the  integral  form  shown 
above t  Applying  the  boundary  conditions,  of  continuity  of  displacements 
and  stresses,  he  arrives  at  six  linear  equations  and  solves  them  for 
functions  f j ,  gl ,  nx ,  f 1 ,  g1,  n1. 

He  then  derives  the  integral  expressions  for  vertical  and  hori¬ 
zontal  displacements  for  PP,  SP,  PS,  SS  and  other  wave  types,  and  then 
solves  them  for  vertical  and  horizontal  displacements  for  P1P2P1, 

P1P2S1  and  P1P2S2.  These  expressions  could  be  evaluated  numerically 


and  compared  with  experimental  data. 


•  ■  -  - 
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II.  3  Travel-Time  Relationships 

Scwhind  et  al ,  (i960)  determine  the  crustal  layering  by  using  PS 
converted  phases  and  the  picture  obtained  thereby,  at  the  sites  men¬ 
tioned  in  their  paper,  favorably  compares  with  that  given  by  other 
methods . 


II.  3.  1  For  an  n-layer  model  the  equation  for  a  wave  travelling 
throughout  as  a  P-wave  as  in  figure  II.  U  is  given  by  Hall  (op,  cit.) 


n-. 


T  =  +  l  ^(^i»vn)(hi  +  hi) 


where  x  is  the  distance  between  shot  point  and  detector 

hj_  is  the  thickness  of  ^th  layer  under  detector 

1 

h-j_  is  the  thickness  of  j_th  layer  under  shot  point 
is  the  compressional  wave  velocity  in  ^th  layer 


and 


n(uj_ , u j  )  =  '  x 
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U: 


1_ 
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.  2  Cos 
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Uj 


ui 


for  all  Uj_  <  Uj 

If  the  conversion  has  taken  place  in  jth  layer  then  the  expression 
Vj  in  n(Vj,Vn)  is  simply  changed  to  Vj ,  indicating  the  conversion  of 
compressional  wave  to  shear  wave  where  Vj  is  the  shear  wave  velocity 
in  jth  layer.  The  same  thing  is  done  if  there  are  more  than  one 
conversion.  This  is  true  for  the  case  for  beds  with  small  dip;  if 
the  dip  is  severe  then  a  dip  factor  enters  into  calculation. 


'f  X  -  >  -;OXC  iC' 
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Mathematically  we  can  see  that 

n(vi,Vj)  >  n(vi>Vj) 

so  that  a  converted  phase  arrives  later  than  the  corresponding  P 
arrival.  The  greater  the  number  of  conversions  the  later  the  event 
arrives.  Physically  it  is  much  more  clear,  since  the  converted  phase 
takes  more  time  than  the  corresponding  P-event  because  of  its  lower 
velocity . 

II.  U  Method  of  Analysis 

The  following  scheme  has  been  adopted  in  the  identification  of 
converted  phases  in  the  present  work. 

1.  All  the  events  were  plotted  on  a  reduced  travel  time  graph 

(T  -  ~)  VS  DISTANCE.  (Fig.  III.  1,  III.  3,  III.  U  and  III,  7) 

0.5 

2.  A  least  square  method,  such  as  given  by  Steinhart  et  al.  (1961), 
was  applied  to  get  velocities  and  intercept  times. 

3.  Based  on  this  velocity  model  along  with  the  well  information  on 
sediments  at  the  shot  sites,  intercept  times  for  various  con¬ 
verted  phases  were  calculated. 

H.  Whenever  these  times  showed  a  reasonable  agreement  with  the 

intercept  times  in  step  3.,  the  event  has  been  identified  as  that 
particular  event. 

5.  The  phase  relationship  of  the  suspected  converted  phase  was  such 
as  to  indicate  a  clear  SV-motion.  This  was  sometimes  not  true. 
(Fig.  III.  8a) 


. 
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60  An  amplitude  measurement  was  made  for  the  later  events  and  compared 
with  compressional  wave  arrival.  (Fig.  III.  9a  and  9b) 

T«  A  qualitative  comparison  of  frequencies  were  also  made.  (Fig. 

II.  5a,  5b,  5c) 

8.  A  few  particle  motion  diagrams  were  constructed  where  the  data 
were  of  high  enough  quality.  (Fig.  III.  8a,  8b,  8c) 

It  can  be  clearly  seen  that  the  time  difference  between  compres¬ 
sional  wave  arrival  and  converted  phase  arrival  gives  the  information 
about  the  thickness  of  the  layer  in  which  conversion  has  taken  place. 
This  is  one  great  advantage  working  with  these  phases.  No  information 
about  the  exact  location  of  shot  point  or  exact  time  of  shot  occurrence 
is  needed  to  derive  the  picture  of  layering.  For  gently  dipping  beds, 
the  dip  could  be  easily  calculated  by  solving  the  different  equations 
involved  for  depths  under  detector  and  shot  point.  Theoretically,  at 
least,  it  is  always  possible  to  find  converted  phases  which  give  in¬ 
formation  about  these  two  quantities  for  a  particular  layer.  In  the 
case  of  deep  curvature  of  the  beds,  as  already  mentioned,  a  dip  factor 
enters  into  calculations. 

II.  U0a,  Some  of  the  criteria  mentioned  in  the  Identification  of 
converted  phases  are  shown  in  a  few  records  reproduced  here,  (Fig, 

II.  5a,  5b,  5c)  In  fig.  II,  5b  are  shown  the  first  arrival  (T^  from 
Moho  and  a  Moho  conversion  (T^),  The  horizontal  component  traces  are 
indicated  by  Radial  and  Transverse  and  the  vertical  component  is 
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Figure  II .  5"b 
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denoted  as  shown .  Radial  motion  can  he  clearly  observed  onPS1  showing 
Sy-motion.  In  fig.  II.  5d  clear  Sy -motion  is  marked  along  with  P-motion, 
whereas  in  fig.  5a  S-arrivals  from  the  subbasement  and  the  intermediate 
layer  respectively  at  62,30  sec.  and  6U.40  sec.  are  shown.  The 
S-arrival  at  69.50  sec,  corresponds  to  Moho  whose  P-arrival  is  also 
shown.  In  all  these  cases  one  can  observe  the  following  characteristics. 

1.  Large  horizontal  motion  (On  S-arrival,  a  large  transverse  motion 
can  be  seen.  Unfortunately  the  radial  trace  was  dead.  The  trans¬ 
verse  motion  with  a  clear  well  marked  spread  velocity  points  to 
Sjj-motion . ) 

2.  Larger  amplitudes ,  and 

3.  Longer  periods. 

The  second  criterion  was  however  not  true  in  all  cases.  But  it 
conforms  with  the  observation  made  by  Andreev  (loc.cit)  that  the 
amplitude  increases  for  converted  phases. 

The  observational  evidence  for  the  generation  of  Sy-wave  by  an 
explosion  is  illustrated  by  a  record  shown  here  in  fig,  II, 5c.  This 
is  a  record  of  the  500  ton  surface  shot  taken  at  a  distance  of  10,000 
feet  from  the  charge  centre.  The  shot  time  is  indicated  by  a  dot  on 
the  left  hand  side  of  the  figure.  The  time  between  two  consecutive 
thick  lines  is  0.20  sec.  The  three  traces  represent  the  radial, 
vertical  and  transverse  motion  respectively.  Approximately  1,20  sec, 
after  the  shot,  the  P-motion  can  be  observed.  Approximately  0,90  sec, 
after  this  motion  the  Sy-type  motion  can  be  observed  on  the  radial 
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trace.  This  represents  the  conversion  of  P  to  S  very  near  the  shot 
point  and  indicates  that  for  this  shot  at  least  it  should  be  possible 
to  observe  events  at  large  distances,  which  have  the  first  part  of  their 
path  as  Sy  type  motion. 
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CHAPTER  III 
EXPERIMENTAL  RESULTS 

III.  1  Records  examined  for  the  purpose  of  the  present  work  belong  to 
two  profiles  as  already  mentioned. 

1.  Bow  City  West.  The  shot  point  was  at  Bow  City*  and  the  recording 
sites  west  of  Bow  City  covering  a  profile  of  about  180  kilometres. 
There  were  fourteen  records  obtained  on  this  profile, 

2.  Suf field  West.  The  shot  point  was  at  Suf field,  and  the  recording 
sites  were  to  the  west,  covering  a  total  distance  of  about  520 
kilometers ,  Only  that  part  of  the  profile  beyond  about  216  kilo¬ 
meters  is  studied  here,  consisting  of  ten  records  from  the  500 
ton  shot.  The  locations  of  these  profiles  are  indicated  on 

fig.  I.  1. 

III.  2  First  Arrivals  From  Bow  City 

Fig.  Ill,  I  is  the  reduced  travel  time  graph  of  various  events 
identified  on  the  Bow  City  records. 

The  spread  velocity  for  first  arrivals  on  this  profile  is  6. l6 
km/sec  for  the  first  seven  records  and  6,50  km/sec  for  the  rest  of  the 
records.  On  two  records  the  first  arrivals  could  not  be  identified 
with  any  confidence  but  immediately  following  is  an  arrival  with  6.50 
km/sec  velocity.  The  two  horizons  represented  by  these  arrivals  are 
respectively  referred  to  as  the  Basement  (6,l6  km/sec)  and  the  Sub- 
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basement  (6.50  km/sec)  with  intercept  times  of  (0.931  *  0.01+5)  sec.  and 
(1.9^7  1  .00U)  sec.  respectively. 

From  the  well  data  the  thickness  of  the  sedimentary  layers  and 
their  velocities  are  known  at  the  Bow  City  Shot  Point.  The  theoreti¬ 
cally  calculated  intercept  time  is  0.870  sec.  Taking  into  account  both 
the  uncertainties  in  the  calculated  and  the  observed  values,  this  is 
considered  to  be  reasonably  good  agreement. 

The  following  picture  emerges  when  the  depth  calculations  are 
carried  out  on  the  basis  of  P-wave  arrivals. 


Vel.  km/sec. 

Depth  km. 

Thickness  km. 

Cretaceous 

3.U1 

1.1+0 

Mississipian 

5.80 

1.1+0 

1,10 

Basement 

6.l6 

2.50 

9.77 

Subbasement 

6.50 

12.27 

23.16 

Intermediate  Layer 

7.25 

35.1+3 

13.27 

Moho 

7.93 

1+8.70 

III.  3  Converted  Phases  from  Basement 

One  can  see  from  the  travel  time  graphs  that  there  are  two  very 
well  correlated  events  with  a  spread  velocity  corresponding  to  6.l6 
km/sec.  Larger  amplitude  and  lower  frequency  of  these  events  point  to 
some  sort  of  an  S-motion  but  the  spread  velocity  is  characteristically 
P- velocity.  This  leads  one  to  conclude  that  the  event  has  been  con¬ 
verted  to  S  at  some  upper  levels.  Intercept  time  considerations  show 
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that  the  conversion  has  taken  place  in  the  sedimentary  layers  as  shown 
in  fig.  III,.  2.  Thus  the  above  two  events  have  been  identified  as 
PPPPS  and  PPPSS.  The  equations  governing  the  two  events  are  given 
below. 


T3  =  7^+  n(v|,v3)h}  +  n(v2,v3)h’  +  n(v1,v3)h1  +  n(v2,v3)h2  (l) 

T5  =  |-  +  n(vj,v3)h'  +  n(V2,V3)h2  +  n(v1,v3)h1  +  n(v1,v3)h2  (2) 


where 


and 


T  is  the  time  taken  by  the  ray  to  travel  from  shot 

point  to  detector  location. 

hj,h2  are  thicknesses  of  the  layers  under  the  shot  point. 

h1,h2  are  thicknesses  of  the  layers  under  the  detectors. 

Vi,V2,V3  are  P-wave  velocities  in  the  three  layers 

v1 ,v2  are  S-wave  velocities  in  the  first  two  layers, 


also 


n(Ui,Uj  ) 


Cos  i. 


'i»J 


th 


i  ui 

where  sin  c  =  —  for  all  U,-  <  U^. 

Uj  } 

It  is  assumed  that  the  following  relationship  holds  between  the  P-wave 
and  S-wave  velocities  in  ^th  layer. 

Vi  =  vi 

The  ray  path  for  P-motion  in  the  above  structure  is  given  by 


T1  =  V7  + 


q(v1,v3)h1  +  n(vj,v3)hj  +  n(v2,v3)h’  +  n(v2,v3)h2 


(3) 


6  .16  km /He 


Fig.  III.  2 
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Fig.  III.  5 
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Solving  (l)  and  (3)  for 
estimates  of  the  values 


h1 ,  one  can  then  solve 
of  hj  and  h2  are  given 


(l)  and  (2)  for  h2>  The 
in  Table  1. 


TABLE  1 


x(km. ) 

T1 ( secs ) 

PPPPP 

T 

12 

sec 

PPPPS 

T2-T! 

.5£C 

hx  (km. ) 

T 

1  3 

PPPSS 

5ec 

T3-T2 

sec 

h2 (km. ) 

176.96 

29.65 

30.22 

0.57 

2.71 

153.59 

25.84 

26.38 

0.54 

2.57 

26.88 

0.50 

2.63 

148.46 

25.04 

— 

— 

— 

25.98 

— 

— 

132.23 

22.39 

— 

— 

— 

23.29 

— 

— 

120.06 

20.42 

20.92 

0.50 

2.38 

21.31 

0.39 

2.05 

110,07 

18.74 

19.28 

0.54 

2.57 

19.68 

0.40 

2.11 

108. l4 

18.46 

18.94 

0.48 

2.29 

19.36 

0.42 

2.21 

100.38 

— 

17.69 

— 

18.12 

— 

86.17 

14.91 

15.41 

0.50 

2.38 

15.61 

0.20 

1.05 

Though  there  is  a  good  amount  of  scattering  in  the  estimation  for 
h}  and  h2 ,  a  general  feature  is  the  thickening  of  both  layers  as  one 
goes  to  the  west  towards  the  mountains.  (See  for  example  Geological 
History  of  Western  Canada  by  the  Alberta  Society  of  Petroleum  Geologists.) 

In  all  the  above  calculations  and  in  the  following  ones  the  lat¬ 
eral  variations  in  velocity  in  the  Cretaceous  formations  has  been 

taken  into  consideration.  It  was  assumed  that  there  is  no  appreciable 

Ml 

change  in  the  velocity  in  Mississippian. 

Ill,  4  Converted  Phases  from  the  Subbasement 

On  the  travel  time  graph  (Fig.  III.  l)  an  arrival  parallel  to 
the  first  arrivals  from  the  subbasement  is  evident.  It  is  not  fully 


. 

. 


understood  at  what  horizon  this  event  originates.  It  seems  to  he  a 
converted  phase  hut  the  times  are  not  consistent  with  a  simple  inter¬ 
pretation. 
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III.  5  Estimate  of  the  Thickness  of  the  Basement 


Since  the  ray  path  travel  time  for  P-motion  from  the  suhhasement 
is  given  hy 


t2  =  ~  +  n(vj,v4)hj  +  n(v1,v4)h1  +  n(v2,v4)h’  +  n(v2>v4)h, 

+  n(v3,v4)(h  +  h’) 


W 


one  can  estimate  the  values  for  (h3  +  h3)  at  different  distances 
w  can  he  written  as 


(h3  +  h3)  = 


T  -  {“  +  n(vj,v4)h1  +  n(v',v4)h2  +  n(v2,v4)h2} 


n(v3>v4) 


(Ua) 


All  the  quantities  on  the  right  hand  side  of  (Ua)  are  known. 

We  can  now  calculate  (h3  +  h3)  at  different  values  of  x.  These  are 
given  in  Table  2. 

The  scatter  in  values  for  (h3  +  h3)  is  not  too  high  compared  to 
the  same  in  h,,  h2,  A  general  feature  of  interest  is  the  thinning 
under  the  mountains  of  the  layer  with  6.1 6  km/sec,  velocity. 
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Table  2 


Bow  City  Arrivals  from  the  Subbasement 
and  their  Interpretation 


X 

K-m 

T—  — - 

Sec 6  •  5 

.19i+h1 

.077h2 

zee. 

T 

Sec 

h3+h3 

^3 

K*ri 

176.96 

1.95 

.526 

mm 

mm 

mm 

mm 

153.59 

1.95 

.U99 

.20 

1 — 1 
OO 

a 

18,08 

8.31 

IU8.U6 

1.95 

— 

— 

— 

— 

— 

132.23 

1.95 

— 

— 

— 

— 

— 

120.06 

1.95 

.U62 

.158 

.893 

19 .81+ 

10.07 

110.07 

1.95 

.^99 

.162 

.81+2 

18.71 

8.91+ 

108. ik 

1.95 

Mk 

.170 

.90 

20.00 

10.23 

100.38 

1.95 

mm 

— 

— 

— 

— 

86.17 

1.95 

.1+62 

.081 

.97 

21,56 

Average  h3 

11.79 

=  9.86 

III.  6  Converted  Phases  from  Intermediate  Layer 

Two  events  were  identified  with  the  same  velocity  as  that  of  the 
intermediate  layer.  The  correlation  from  record  to  record  is  good. 

On  the  basis  of  intercept  times  these  two  events 9  in  the  order  of 
occurrence,  have  been  respectively  identified  as 

1.  conversion  in  basement  and  sediments 

2.  conversion  in  basement  and  sediments,  and  conversion  in  subbase- 
ment  only  under^shot  point. 

Table  3  gives  the  estimations  of  thicknesses  of  various  layers 

jLXX 

based  on  these  conversions  from  Intermediate  Layer. 
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Table  3 


Intermediate  Layer  Arrivals 


X 

kyr\ 

Ti 

ppppp 

Sec. 

T 

12 

SPPPS 

Sc-C. 

T2-T! 

h3+h3 

/Cvn 

£3 

T 

1  3 

SSPPS 

T  -T 

3  2 

See. 

h4 

M  yr\ 

153.59 

27.18 

— 

mm 

... 

mm 

mm 

ll+8. 1+6 

26.1+9 

30.89 

1+.1+0 

19.10 

9.33 

35.29 

1+ .  1+0 

26.39 

132.23 

21+.21+ 

28.61+ 

1+  .1+0 

19.10 

9.33 

32,99 

1+  .35 

26,06 

120.06 

22.57 

26.92 

1+.35 

18.76 

8.99 

31.31 

h.39 

26.32 

110.07 

21.18 

25.53 

1+.35 

18.76 

8.99 

29.85 

1+.32 

25.87 

.100.38 

19.81+ 

21+.19 

1+.35 

18.76 

8.99 

28.59 

1+  0 1+0 

26.39 

Aver.  = 

9.13 

Aver.  = 

26.20 

The  following  picture  emerges  from  these  calculations. 

Thickness 

Velocity 

Shot  Point 

km 

km/sec , 

*Cretaceous 

1.1+ 

3.1+7 

*Mississippian 

1.1 

5.80 

Basement 

9.77 

6 . 16 

Subbasement 

26.20 

6.50 

Detectors 

Cretaceous 

2.1+8 

1+.01 

Mississippian 

ro 

0 

0 

0 

5.80 

Basement 

9.13 

6,16 

Subbasement 

20.95 

6,50 

There  is  an  apparent 

discrepancy  in  the 

value  for  the  subbase 

*  From  the  well  data 
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ment  thickness  under  the  shot  point  as  calculated  in  two  different  ways 
i.e.  from  first  arrivals  and  from  conversions  from  the  intermediate 
layer 0  This  is  mainly  due  to  the  small  quantities  involved  in  the 
calculations  for  these  converted  phases.  The  times  are  small  and  the 
n's  are  small.  A  change  even  in  third  decimal  place  in  n’s  causes  a 
change  of  3  -  ^  km  in  thickness  estimates. 

III.  T  First  Arrival  Data  From  Suffield  500  Ton  Profile 

In  fig.  III.  3  and  III.  4  are  given  the  reduced  travel-time 
graphs  of  some  events  observed  on  500  ton  records. 

The  first  arrivals  from  500  ton  profile  have  a  spread  velocity 
of  about  8.37  km/sec.  The  correlation  from  record  to  record  is  fairly 
satisfactory.  The  intercept  time  is  9*25  *  0.l4  sec.  Following  this 
event  a  number  of  other  events  were  identified  and  plotted  on  the 
graph.  The  event  with  velocity  7.11  km/sec  corresponding  to  the 
intermediate  layer  has  an  intercept  time  of  5.60  ±  0.12  sec.  The 
following  table  gives  the  intercept  times  and  velocities  for  both 

Table  4 


Velocity  Km/Sec.  Intercept  Time  Sec, 


Moho 

Inter  Layer 

Moho 

Inter 

Layer 

Bow  City 
Profile 

7.93  * 

.05 

7.24  *  .002 

8,64  *  .10 

5.99 

*  ,05 

500  Ton 

Profile 

8.37  ± 

.04 

7.11  *  .004 

9.25  *0.14 

5.60 

*0,12 
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profiles  for  these  two  horizons. 

The  correlation  for  the  intermediate  layer  for  the  500  ton  pro¬ 
file  is  not  good,  with  the  result  that  the  event  has  been  shown  on  the 
plot  by  two  line  segments.  One  of  the  segments  fits  in  with  the 
remainder  of  the  Suf field  West  profile.  These  two  lines  have  differ¬ 
ent  velocities  and  intercept  times,  and  if  the  correlation  is  correct, 
they  indicate  strong  dip  on  the  top  of  the  7«1  km/sec.  layer. 

There  is  no  arrival  from  the  basement.  However,  there  are  a 
number  of  arrivals  from  the  subbasement.  It  was  therefore  assumed  that 
the  basement  layer  is  either  missing  or  very  thin  under  the  Suffield 
shot  point,  a  circumstance  to  be  confirmed  later  from  the  converted 
phases . 

The  thicknesses  of  sediments  and  velocities  in  these  formations 
under  the  Suffield  shot  point  are  known  with  fair  accuracy  from  well 
data. 

The  thickness  of  the  6.50  layer  was  calculated  on  the  assumption 
that  the  basement  layer  is  absent  at  Suffield.  The  value  is  31.50  km. 
Similarly  the  depth  to  the  intermediate  layer  was  calculated.  The 
following  picture  emerges  from  these  calculations. 


Under  Suffield  Shot  Point 


Thickness 

km 


Velocity 
km/sec . 


^Cretaceous 


1.30 


2,11 


*  From  well  velocity  survey 
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Thickness  Velocity 

km  km/sec. 

*Mississippian  0.84  5.80 

Subbasement  31.50  6.50 

Intermediate  Layer  10.00  7.11 

Depth  to  Moho  43.64  km. 


Velocity  below  Moho  8.34  km/sec. 
III.  8  Converted  Phases  from  the  Mohorovicid 


The  first  of  these  phases  occurs  just  about  2.30  seconds  after 
the  first  arrival.  This  has  been  identified  on  the  basis  of  intercept 
times  as  the  conversion  taken  place  at  the  basement-subbasement  bound¬ 
ary,  and  is  shown  in  Fig.  III.  5»  The  time  difference  between  the 
first  arrival  which  is  P-motion,  and  this  phase,  therefore  gives  the 
thickness  of  the  formations  overlying  the  subbasement.  The  following 
relations  hold  for  these  two  phases. 


T1  = 


V, 


n(v1,v6)(h1  +  hj)  +  n(v2,v6)(h2  +  h’) 


n(v3,v6)(h3 


+  n(v4,v6)(h4  +  h4)  +n(v5,  v6)(h5  +  h5) 

T2  =  ~  +  n(v1,V6)(h1  +  h’)  +  n(v2,V6)(h2  +  h2)  +  ri(v3,V6)(h3  +  h3) 
V6 

+  ^\»V6)(h4  +  O  +  r'(V5*V6  )(h5  +  h5) 


*  From  well  velocity  survey 
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where  V1 ,  .  .  .  Vg  are  P-wave  velocities  in  the  layers  and  v}  ,  v2  ,  v3 
are  S-wave  velocities  in  the  first  three  layers .  Clearly  the  differ¬ 
ence  (T2  -  T  )  =  AT  gives  the  time  taken  by  the  wave  to  travel  in  the 
first  three  layers.  The  thickness  of  the  sediments  is  known  under  the 
Suffield  shot  point  and  under  the  detectors.  If  is  the  combined 
thickness  of  the  three  layers  under  detectors  and  hJ  the  same  under 
shot  point,  the  calculations  yield  as  shown  in  table  6,  a  figure  of 
average  thickness  16.U8  km  for  (H1  +  h|)„  Two  or  three  of  the  detector 
positions  for  500  ton  profile  lie  very  close  to  those  for  Bow  City 
profile.  The  distance  between  Bow  City  and  Suffield  is  126.50  km. 

Some  detector  positions  for  both  profiles  are  given  below  in  table  5<= 

Table  5 

Suffield  Bow  City  x  +  126.50 

x 

Kvry  A'w,  ™ 

216.15  86.17  212.67 

276.85  1U8.U6  27^o96 

28U.71  153.59  280.09 

298.62  176.96  303.U6 


From  the  above  table  it  can  be  easily  seen  that  at  least  two  of  the 
detector  positions  are  within  2-3  km.  for  both  profiles;  two  of  the 
positions  are  at  least  within  5  km.  for  both  profiles.  The  informa. 
tion  at  these  places  is  already  known  from  Bow  City  profile. 

From  the  Bow  City  profile,  therefore,  we  take  the  following 
figures  for  the  thickness  of  sediments  and  6,l6  layer. 


■ 


dSi  si  Jb.ie..  r":i  ■  b  y  riO  w'<€  =»■>' 


...  .  .  • 


Sediments 


U . 30  km 


3b 


Basement 


9.86  km 


Thickness  of  sediments  under  Suf field  2,1*1  km 


Total  (average  figure) 


l6o30  km 


The  corresponding  figure  as  calculated  from  the  converted  phase 
for  500  ton  profile  is  l6 0 U8  km.  This  is  evidence  for  extreme 
thinness  (or  absence)  of  the  6.l6  layer  under  Suf field,  as  shown  in 
Fig.  III.  5.  This  was  already  pointed  out  in  the  interpretation  under 
first  arrival  data. 

The  next  converted  phase  occurs  approximately  8  secs,  after  the 
first  converted  phase  and  approximately  10  secs,  after  the  first  arrival 
from  the  Moho.  On  the  basis  of  intercept  times  this  has  been  inter¬ 
preted  as  the  conversion  taken  place  in  the  layers  up  to  and  including 
the  subbasement.  This  therefore  gives  information  about  the  subbasement. 

The  following  relation  holds  for  it. 


+  n(v4,Y6)(h4  +  h4)  +  n(V5,V6)(h5  +  hg) 


where  v4  is  the  S-wave  velocity  in  the  subbasement. 
Clearly  (T  -T2)  =  AT  is  given  by 


at  =  [n(v4#v6)  -  n(v4,v6)]  (h4  +  h4) 


or 


/,  I  n  _ AT _  _  AT 

4  4  ~  n(v  ,v  )  -  n(v  ,V")  An 

4  6  4  6 


, 


.  •  .  7'  i  r1*  ^  ' 
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Since  all  the  quantities  on  right  hand  side  are  known  one  can 
calculate  (h4  +h4 ) .  These  estimates  are  given  in  table  6.  The  aver¬ 
age  value  for  (h4  +  h4)  is  thus  51.62  km.  But  thickness  of  the  sub- 

l 

basement  under  the  Suffield  shot  point  is  31c 50  km.  as  calculated 
previously.  This  leaves  a  thickness  of  20.12  km.  for  the  subbasement 
under  the  detectors.  It  should  be  remembered  here  that  the  detectors 
are  located  quite  far  away  from  Suffield.  Some  detector  positions 
are  in  the  foothills  while  others  are  located  deep  in  the  mountains. 
This  then  says  that  the  subbasement  is  getting  thinner  under  the 
mountains.  Further,  this  figure  favorably  compares  with  the  thickness 
calculated  for  the  same  horizon  in  the  Bow  City  profile  in  III.  6,  the 
figure  being  20.95  km. 

The  next  Converted  phase  ocCubs  approximately  i+.50  secs,  after 
the  second  conversion  and  lU.50  secs,  after  the  first  arrivals  from  the 
Moho.  This  has  been  interpreted  again  on  the  basis  of  intercept  times, 
as  the  conversion  taken  place  in  all  layers  above  Moho.  The  equation 
for  this  is  given  below. 

Tg  -  +  n(v1,V6)(h1  +  hj)  +  n(v2,V6)  (h2  +  h2)  +  n(v3,v6)(h3  +h3) 

v  6 

+  n(v4,V6)(h4  +  h’)  +  n(v5,V6)(h5  +  hj) 

The  difference  AT  =  (Tg  -  T7 )  is  thus 


AT  =  T9  -  T7  =  [n(v5,V6)  -  n ( V5 ,V6 ) ]  (h5  +  hg) 


■  r 
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or  (h5  +  h’)  =  ~ 

5  b  z\n 

where  are  S-  and  P-wave  velocities  in  Intermediate  Layer  and 

» 

h59  h5  are  the  thicknesses  of  this  layer  under  detectors  and  shot 
point  respectively. 

This  thus  gives  the  thickness  of  the  intermediate  layer  under  the 
shot  point  and  detectors  combined.  The  average  figure  for  this  is 
25.23  km.  But  the  thickness  of  this  layer  under  shot  point  as  given 
under  the  interpretation  of  first  arrival  data  is  about  10  km.  This 
leaves  a  thickness  of  about  15  km.  under  the  detectors.  This  then 
says  the  intermediate  layer  is  thicker  under  mountains. 

The  following  picture  for  crustal  layers  thus  arises  from  con¬ 
verted  phases  under  the  detectors. 


Sediments 

U.30 

km 

Basement 

9.86 

km 

Subbasement 

20.12 

km 

Intermediate  Layer 

15  °23 

km 

Depth  to  Moho 

49  o71 

km 

III.  9  Converted  Phases  from  Intermediate  Layer 

The  first  of  these  occurs  round  about  8  secs,  after  the  first 

event  from  the  intermediate  layer.  On  the  basis  of  intercept  times 

a. 

this  has  been  interpreted  as ^conversion  taken  place  above  the  sub¬ 
basement  only  under  the  shot  point  as  shown  in  fig.  III.  6.  The 
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rest  of  the  path  it  travels  as  P-wave,  The  equation  for  this  is  given 
below  as  T&,  T.^  is  the  equation  for  the  wave  having  the  same  path, 
except  that  it  travels  as  P-wave  throughout, 

t3  =  n(v1,v5)(h1  +  hj)  +  n(v2,v5)(h2  +  h2)  +  n(v3,v5)(h3  +  h3) 

+  n(v4,v5)(h4  +  h*) 

T6  =  +  n(visv5)h  +  n(v1,v5)hj  +  n(v2,v5)h2  +  n(v2,v5)h2 

5 

+  n(v3,v5)h’  +  n(v3,v5)h3  +  n(v4,v5)h4  +  n(v48v5)h4 

4 

at  =  t6  -  t3  =  z  [n(vi#v5)  -  n(v.sv5)]h! 

h^  is  the  thickness  under  shot  point  of  ^th  layer.  Since  the  thick¬ 
nesses  of  and  velocities  in  sediments  under  Suffield  are  known  one  can 
calculate  the  thickness  of  subbasement  which  gives  an  average  value 
of  30,26  km.  This  agrees  favorably  with  the  figure  we  obtained  under 
the  interpretation  of  first  arrival  which  was  31<>50  km. 

Following  this  there  is  another  arrival  with  the  same  spread 
velocity  occurring  about  4  secs,  after  the  above  arrival  and  about  12 
secs,  after  the  first  event  from  the  intermediate  layer.  It  has  been 
identified  again  on  the  basis  of  intercept  times  as  a  conversion  taking 
place  in  the  sediments  and  basement  on  one  path  and  a  conversion  in  the 
subbasement  only,  on  the  other  part  of  the  path.  The  equation  for  this 


is  given  by 


:  :/ 
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t8  =  +  n(v1,V5)(h1  +  h] )  +  n(v2SV5)(h7  +  h2)  +  n(v39V5)(h3  +  h3) 

+  n(v4,V5)h’  +  n(V49v5)h4 
3 

AT  =  Te  -  T6  =  2  t n (vi 9V5 )  -  n(v.9v5)]h.; 

i=l 

This  thus  gives  a  relationship  between  the  time  under  and  the  thick¬ 
nesses  of,  the  sediments  and  basement  under  the  detectors.  Table  7 
gives  the  thicknesses  calculated  for  various  layers  on  the  basis  of  the 
converted  phases  from  the  intermediate  layer.  Knowing  the  thickness  of 
sediments  one  can  then  calculate  the  thickness  of  basement.  The  thick¬ 
nesses  thus  calculated  at  different  distances  are  8.49,  8 .44,  8.49 
with  an  average  cf  8.47  km.  The  average  thickness  for  the  same  layer 
under  Bow  City  profile  was  9*86  km. 

The  crustal  picture  down  to  the  intermediate  layer  from  these 
calculations  is  as  follows. 
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III.  10  S-Phases  and  Lg  Phase 


In  addition  to  the  various  converted  phases  discussed  above, 
S-phases  from  three  horizons,  the  Subbasement,  Intermediate  Layer  and 
Mohotovicic ,  were  identified.  Fig.  Ill,  7  is  a  travel-time  graph  for 
these  arrivals.  Table  8  gives  P  and  S  velocities  and  Poisson’s  ratio 
for  these  layers. 


Table  8 


Horizon 

Vp 

Mi  t<rr  -  - 

Vs 

A 

Vp/Vs 

a 

Subbasement 

6.50 

3.75 

1.73 

.21+9 

Intermediate  Layer 

7.11 

1+.53 

1.57 

.159 

Moho 

8.37 

U.73 

1,76 

.263 

following  table 

gives  the  distances  and  time  for^Lg  Phase, 

Table  9 

£td\.l)'or\ 

Distance  (km,) 

Time  (secs,) 

Velocity 

km/sec . 

Kilo 

218.15 

61,12 

3.57 

Tango 

21+3.55 

67.65 

3.60 

Juliet 

251+.1+5 

— 

- 

Hotel 

270.62 

72.21 

3.75 

Sierra 

275. 3U 

76.38 

3.60 

Lima 

283.21 

77.15 

3.67 

Papa 

300,11 

83.35 

3.60 
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This  phase  is  particularly  well  displayed  on  the  record  at  a  distance 
of  218  km.  where  a  clear  onset  and  a  dispersive  wave  train  can  he  seen. 
The  presence  of  this  phase  is  another  clear  indication  that  S  waves 
are  generated  by  an  explosive  source  at  least  in  this  particular  case. 


III.  11  In  table  10  are  enumerated  the  calculated  and  observed 
intercept  times  for  all  the  events  observed  on  both  profiles. 

Table  10 
Intercept  Times 


Event 

a.  Suffield  500  Ton  Profile 

1.  Mohorovicid 

Ti 

T2 

T7 

2.  Intermediate  Layer 

T3 

t6 

t8 


Calculated  Observed 

Time  Time 

Sec  sec. 


9.2b 

9.25 

± 

,1b 

11.98 

11,10 

± 

.33 

19.31 

18.75 

± 

1.35 

22.77 

23.25 

± 

0.55 

5.63 

5.60 

± 

.12 

13. b2 

15.5b 

± 

1,10 

18,23 

18,92 

b,  Bow  City  Profile 


T, 

.87 

.93 

± 

.05 

T2 

1.97 

1.95 

± 

t3 

l.b3 

1.32 

± 

.07 

t5 

2,b2 

2,b6 

± 

.07 

T6 

5.89 

5.99 

± 

.05 

T7 

8.5b 

8,6b 

± 

.10 

T 

18 

T9 

10. bO 

10.17 

± 

,08 

13.95 

lb.  50 

± 

.319 
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III.  12  Particle  Motion  Diagrams 


In  several  cases  where  the  event  was  of  high  quality  with  a  good 
horizontal  motion  recorded,  the  particle  motion  diagram  was  drawn. 

These  diagrams  are  shown  in  Figs.  III.  8a  and  III.  8b,  Polarisation 
of  PS  phase  is  not  particularly  clear  from  these  diagrams,  even  though 
there  is  a  large  horizontal  component  of  motion.  The  only  bright 
feature  about  these  diagrams  for  PS  phases  is  that  they  indicate  a 
predominantly  horizontal  motion  with  a  phase  difference  between  ver¬ 
tical  and  horizontal  components.  The  following  factors  have  to  be 
borne  in  mind  while  collecting  the  records,  in  order  to  draw  useful 
particle  motion  diagrams. 

1.  Orientation  of  horizontal  geophones, 

2,  All  the  other  factors  mentioned  under  amplitude  measurements  below, 

III.  13  The  amplitudes  of  some  of  the  events  were  measured.  The 
ratios  of  amplitudes  of  the  converted  phases  from  the  Mohorovicid  to 
that  of  the  P-event  are  plotted  against  the  distance  in  Fig,  9a  and 
the  same  for  the  basement  and  subbasement  in  Fig.  III.  9b,  A  large 
scatter  in  the  data  is  evident  and  it  did  not  seem  worthwhile  to  make 
detailed  comparisons  with  theoretical  values.  If  more  data  were 
available  it  might  be  possible  to  obtain  some  smoothed  estimates  of 
the  amplitude  ratios  and  hence  obtain  results  which  could  be  compared 
in  a  meaningful  way  to  calculated  values. 

All  the  amplitude  measurements  are  in  arbitrary  units  and  these 
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measurements  are  not  quantitative  by  any  means ,  No  comparison  could 
be  made  as  there  were  few  measurements  as  mentioned  already  and  no 
conclusions  drawn.  Certain  features  which  are  to  be  taken  care  of 
while  the  seismic  surveys  are  done  are  the  following, 

1,  Information  about  surface  relief  where  the  geophone  is  located  is 
necessary. 

2.  Proper  settings  for  attenuation  and  filter  settings  to  be  noted 
down  clearly  on  the  records . 

3o  ^void  as  far  as  possible  the  power  lines. 

4.  Preferably  all  records  and  all  traces  should  be  recorded  on  tape. 
Once  the  records  are  on  tapes  it  is  easier  to  apply  crosscorrelation 
techniques  and  filtering  techniques  to  reduce  the  noise  and  get  a  more 
accurate  amplitude  measurement. 
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CHAPTER  IV 
CONCLUSIONS 


An  attempt  thus  has  been  made  to  interpret  various  later  events 
observed  on  two  profiles.  The  calculations  carried  out  based  on  the 
formulae  given  by  Hall  (loc.  cit.)  are  in  some  cases  in  reasonably 
good  agreement,  while  in  some  cases  not,  but  the  overall  picture  of 
crustal  layering  interpreted  by  this  method  is  fairly  convincing. 

The  depth  to  the  Moho  under  Suffield  is  about  kk  km. , under  Bow  City 
4o  km.  and  under ^mountains  49  km.  Remembering  that  Bofe/  City  is  west 
of  Suffield,  this  picture  of  a  deeper  Moho  towards  the  west  is  in 
agreement  with  the  conception  of  thickening  crust  under  the  Rockies, 
although  this  thickening  is  not  as  great  enough  to  suggest  a  root 
under  the  Rockies.  Further  this  figure  of  49  km.  must  have  rather 
large  confidence  limits  assigned  to  it  and  hence  the  value  is  con¬ 
sidered  to  be  in  reasonable  agreement  with  other  estimates  in  the 
same  area. 

Taking  layer  by  layer  we  see  first  that  the  sediments  are 
getting  thicker  near  the  foothills.  (Figs.  IV.  1,  IV.  2)  This 
agrees  with  the  geological  evidence  of  the  occurrence  of  Mississip- 
pian  outcrops  in  the  mountains.  Secondly  the  basement  itself  is 
getting  thinner  under  the  mountains.  Thirdly  the  subbasement  does 
the  same  and  lastly  the  intermediate  layer  gets  thicker  under  mountains. 
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Certain  discrepancies,  for  example  in  the  thickness  of  the  suhhasement 
under  the  Bow  City  shot  point,  did  occur  in  these  calculations ,  This 
is  because  of  the  very  small  quantities  involved  in  many  of  the  calcu¬ 
lations,  The  functions  n(U^,Uj)  are  very  small  and  the  difference  An 
is  still  smaller  with  the  result  a  change  even  in  the  third  decimal 
place  affects  the  thickness  by  at  least  3-^-  km.  This  ultra  sensitive¬ 
ness  of  these  functions  makes  it  necessary  to  know  the  true  velocity  in 
a  given  layer  to  a  high  degree  of  accuracy.  This  was  not  achieved  in 
this  case  due  to  lack  of  reverse  control  on  the  various  horizons. 

This  has  been  a  big  handicap  in  these  investigations. 

Comparison  of  the  theoretically  calculated  intercept  times  with 
the  observed  ones  shows  that  for  the  first  arrivals  these  two  agree 
reasonably  well.  For  example  the  first  arrivals  on  the  Suffield  500 
ton  profile  have  a  calculated  intercept  time  of  9«2U  sec,  against  the 
observed  9.25  sec.  For  Bow  City  the  basement  arrivals  have  a  calculated 
intercept  time  of  0,87  against  the  observed  0,88  (taking  the  lower 
limit,)  For  the  subbasement  the  figure  is  1,97  against  an  observed 
1.95.  The  same  figure  for  intermediate  layer  for  Suffield  500  ton  is 
5,63  against  the  observed  5.60,  Most  of  these  calculated  times  thus 
are  either  reasonably  close  to,  or  within  the  90%  confidence  limits  of, 
the  observed  ones.  Mention  has  to  be  made  here  that  all  the  intercept 
times  have  been  calculated  on  the  basis  of  first  arrival  data.  Hence 
any  discrepancy  in  the  calculation  of  intercept  times  for  later  arrivals 
can  be  easily  understood  as  a  consequence  of  applying  the  results  of 
one  event  to  another  which  are  completely  different  in  character. 
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